import os, sys
from pylab import *
from scipy.ndimage import gaussian_filter1d
from PhotozSN_Disp import *

#sn_sample = PhotozSNeDisp(  n_mock=100,\
#                            sample_size=350000,\
#                            mock_sn_dir="SNR15_Nfilt3_main",\
#                            sn_dist_outroot="SNR15_Nfilt3_main")

#sn_sample.estimate_disp()

#sn_sample = PhotozSNeDisp(  n_mock=100,\
#                            sample_size=35000,\
#                            mock_sn_dir="SNR15_Nfilt3_main_35000",\
#                            sn_dist_outroot="SNR15_Nfilt3_main_35000")

#sn_sample.estimate_disp()

#sn_sample = PhotozSNeDisp(  n_mock=100,\
#                            sample_size=12000,\
#                            mock_sn_dir="SNR15_Nfilt3_deep",\
#                            sn_dist_outroot="SNR15_Nfilt3_deep_12000")

#sn_sample.estimate_disp()


#sn_sample = PhotozSNeDisp(  n_mock=100,\
#                            sample_size=100000,\
#                            mock_sn_dir="SNR15_Nfilt3_main_100000",\
#                            sn_dist_outroot="SNR15_Nfilt3_main_100000")

sn_sample = PhotozSNeDisp(  n_mock=100,\
                            sample_size=400000,\
                            mock_sn_dir="SNR15_Nfilt3_main_400000",\
                            sn_dist_outroot="SNR15_Nfilt3_main_400000")

sn_sample.estimate_disp()
